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ABSTRACT 

We present a quantitative morphological analysis using Hubble Space Telescope (HST) 
NICMOS Hiqq- and ACS /775-band imaging of 25 spectroscopically confirmed sub- 
millimetre galaxies (SMGs) which have redshifts between z = 0.7-3.4 (z = 2.1). Our 
analysis also employs a comparison sample of more typical star-forming galaxies at 
similar redshifts (such as Lyman Break Galaxies) which have lower far-infrared lu- 
minosities. This is the first large-scale study of the morphologies of SMGs in the 
near-infrared at ~ 0.1" resolution (< lkpc). We find that the half light radii of the 
SMGs (r^=2.3=b0.3 and 2.8±0.4kpc in the observed I- and i^-bands respectively) and 
asymmetries are not statistically distinct from the comparison sample of star-forming 
galaxies. However, we demonstrate that the SMG morphologies differ more between 
the rest-frame UV and optical-bands than typical star-forming galaxies and interpret 
this as evidence for structured dust obscuration. We show that the composite ob- 
served i7-band light profile of SMGs is better fit with a high Sersic index (n ~ 2) 
than with an exponential disk suggesting the stellar structure of SMGs is best de- 
scribed by a spheroid/elliptical galaxy light distribution. We also compare the sizes 
and stellar masses of SMGs to local and high-redshift populations, and find that the 
SMGs have stellar densities which are comparable (or slightly larger) than local early- 
type galaxies, but comparable to luminous, red and dense galaxies at z ~1.5 which 
have been proposed as direct SMG descendants, although the SMG stellar masses 
and sizes are systematically larger. Overall, our results suggest that the physical pro- 
cesses occuring within the galaxies are too complex to be simply characterised by the 
rest-frame UV/optical morphologies which appear to be essentially decoupled from all 
other observables, such as bolometric luminosity, stellar or dynamical mass. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift - sub- 
millimetre 



1 INTRODUCTION 

Around 60% of the stellar mass in the local Universe is con- 
tained within early- type and elliptical galaxies, which sit 
on a tight "red sequence" in the colour magnitude diagram 
(ISandage fe Visvanathanll 19781 : iBower et al.lll992l : iBell et al.1 
2003). Early-type and elliptical galaxies follow well known 
scaling relations (the fundamental plane) and exhibit sys- 
tematic correlations between the absorption line strengths 



and velocity dispersion, (a). This age-cr relation is such that 
the most massive (a ~400kms -1 ) galaxies appear to have 
formed their stars ~ 10-13 Gyr ago. By contrast the mean 
age of the lower dispersion galaxi es (a ^50 km s -1 ) formed 
more recently; ~4Gyr ago (e.g. ISmith et al.l l2007h . This 
suggests that the stars in giant re d galaxies formed early 
in the history of the Universe (eg. Nela n et al.l [2005). Us- 
ing deep near-infrared imaging, it has also become possi- 
ble to extend the selection of red galaxies to higher red- 
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shift and identify massive, relatively old galaxies at z ~ 1.5 
which could be cons idered the progenitors of the local ellip- 
tical population fee. Ivan Dokkum et alJliooS : ICimatti et al.l 
2008). Clearly to probe this evolutionary sequence further, 
direct observations of the formation of the most massive 
galaxies at high-redshift are required. However, this has 
turned out to be a non-trivial exercise since the most ac- 
tively star- forming, massive galaxies at z >2 are also the 
most dust obscured ([Dole et al.l 1200^1 ; IPapovich et al. 2004; 
iLeFloc'h et"alll2009h . 

Nevertheless, mid- and far- infrared surveys (particu- 
larly those made with the 850/mi SCUBA camera on the 
JCMT and more recently with the Spitzer Space Telescope 
at 24/im) have begun to resolve the most highly obscured 
populations into their constituent galaxies, determine their 
contribution to the energy density in the extra-galactic far- 
infrared/sub- mm back ground, and chart t hat history of mas- 
sive galaxy formation ([Smail et al.l 12002k ICowie et al.l 12002; 
iLe Floc'h et~alll2005h . Extensive, mult i- wavelength follow- 
up has shown that these heavily dust-obscured, gas-rich 
galaxies lie predominantly at high redshift (z ~ 2; e.g. 
IChapman et al.l l2003aL l2005h , with bolometric luminosities 
of ^>1O 12 L and star-formation rates of order 700 M yr _1 . 
It has therefore been speculated that SM Gs are the progen- 



itors of luminous elliptical galaxies (e.g. iLillv et al 



Genzel et al.ll2003l: iBlain et alJliool ISwinbank et al 



1999 



2006 



Tacconi et al.l 12008). 



With the redshift distributions and contributions to the 
cosmic energy density of ultra- luminous galaxies reasonably 
constrained, the next step is to study the evolutionary his- 
tory of SMGs, and to determine how they relate to lower 
luminosity galaxies. Indeed, given the apparently rapid evo- 
lution in the space density of ULIRGs from z ~ 2 to z — 
([Chapman et alJliooEl : ILe Floc'h et al.ll2005h , one key issue 
is to understand the physical processes which trigger these 
far-infrared luminous events. Indeed, the mechanism respon- 
sible for these vigorous starbursts is still uncertain. Analogy 
to local ULIRGs would argue for merging as the trigger, al- 
though secular bursts in massive gas disks is also conceivable 
and indeed recent theoretical interest has stressed the im- 
portance of cold flows in high-redshift star formation (e.g. 
iGenel et afl l2008h . The suggestion that SMGs have com- 
pact disk-like ga s reservoirs (Rj / 2 < 2 kpc) with "maxi- 
mal starbursts" (|Tacconi et al.ll2008h hints that SMGs are 
scaled-up versions of the local ultra-luminous galaxy pop- 
ulation, which are u sually associated with merger activity 
(jTacconi et alJl2002h . 

In order to test the connection between SMGs, lower- 
luminosity star- forming galaxies at high redshift, as well 
as local ULIRGs, we have obtained high resolution Hub- 
ble Space Telescope (HST) imaging of a sample of spec- 
troscopically confirmed SMGs at z =0.7-3.4. By neces- 
sity, most morphological studies of high redshift galax- 
ies to date have been performed at optical wavelengths 
which probe the rest-frame UV (IChapman et al. 2003b; 
Web b et al.ll2003l: IConselice et al.l [2003b. 20081: ISmail et al.l 



2004; Law et al . 2007b). However the rest-frame UV is dom- 



inated by radiation which traces the brightest, active star- 
forming regions rather than the bulk of the stellar pop- 
ulation and can lead to late-type galaxies being classi- 
fied as irregular systems ([Dickinson et al.l l2000l ; iThompsonl 
2003; Goldader et al.|[2002i ). Nevertheless, in a recent study, 



lLaw et all |2007b) conducted a detailed analysis of the rest- 
frame UV morphologies of a large sample of UV/optically 
selected star- forming galaxies at z ~l-3 in the Hubble Deep 
Field North (HDFN) and find evidence that dusty galax- 
ies have more nebulous UV morphologies than more typ- 
ical sources, but otherwise conclude that UV morphology 
is statistically decoupled from the majority of physical ob- 
servables (such as stellar or dynamical mass, gas fraction 
or star- format ion rate). Here we aim to extend this work 
to include the rest-frame optical emission to test whether 
there are key differences in the morphologies at longer wave- 
lengt hs (as suggested by imaging of low redshift ULIRGs; 
e.g. iGoldader et al.ll2002h . We have therefore assembled a 
sample of 25 SMGs with both HST ACS /-band and NIC- 
MO S if -band observations. We determine the basic morpho- 
logical parameters of this sample of SMGs, as well as search 
for signs of tidal features and major mergers. To baseline 
our analysis we also use a sample of 228 optically selected 
star-forming galaxies at z ~ 2-3 (of which 53 have also been 

observed in if -band with HST). 

We use a WMAP cosmology (Soe rgel et aDl2Q03h with 
ft A =0.73, On=0.27, and H =72kms _1 Mpc -1 . In this cos- 
mology, at z =2.1 (the median redshift of our SMG sample), 
0.1" (the typical resolution of our observations) corresponds 
to a physical scale of 0.8 kpc. All quoted magnitudes are on 
the AB system unless otherwise noted. 



2 OBSERVATIONS, REDUCTION AND 
ANALYSIS 

2.1 SMG Sample 

We used HST NICMOS during Cycle 12 (PID: 9506) to ob- 
serve 23 galaxies spanning the r edshift range 0.7 < z < 3.4 
from the spectroscopic catalog of lChapman et all ([2005). We 
also include in our analy s is two spectrosc opically confirmed 



SMGs from iBorvs et al.l d200^) fsee also IPope et al.1 12005) 



Chapman et al . (2005) sample and 



which are part of the 
which lie in the GOODS-North (GOODS-N) field and were 
observed with NICMOS during program PID: 11082 (Con- 
selice et al. 2010 in prep). Therefore the final sample consists 
of 25 SMGs with secure spectroscopic redshifts, and local- 
ized via VLA radio imaging. 

The majority of the NICMOS targets also have optical 
coverage from HST as part of the same program from Cy- 
cle 12, although we also include archival data in the ELAIS 
field (PID: 9761) and GOODS Legacy program. Briefly, 20 
galaxies have HST ACS imaging, four galaxies have WFPC2 
observations, and one galaxy has STIS imaging. Unless oth- 
erwise stated, we will refer to the ACS, WFPC2, and STIS 
images as "optical" , and specifically the /-band filter. In the 
case of the lone STIS image, no filter was used (hence the 
wavelength range covered is ^5400-1/im). We note that the 
sample spans a redshift range z=0.7-3.4 (median redshift 
z =2.1), 850/xm flux (5*850 =3-15 mJy) and /-band magni- 
tude range (i=21.5-26.5) which is representative o f the par- 
ent sample of 73 SMGs in IChapman et al.l ([20051 ). We also 
note that there are roughly equal numbers for each Az = 1 
bin (Table 1). 

The near-infrared observations were made using the 
NIC2 camera and the F160W filter. We employed the stan- 
dard spiral dither pattern under LOWSKY conditions. Each 
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Figure 1. High resolution HST optical (ACS/WFPC2/STIS) and near-infrared NICMOS imaging of SMGs. For each galaxy, the left 
hand panel denotes the optical image whilst the right hand panel denotes the 1.6/mi (if —band) image. Each thumbnail is scaled such 
that they each cover 40 kpc at the redshift of the SMG. The small cross denotes the location of the centroid of the radio emission. 
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Figure 2. True colour HST IH-b&nd images of the SMGs in our sample showing the range of colours and morphological mix within the 
sample. Each image is 40 kpc at the redshift of the galaxy. 





Figur e 3. C omparison true colour HST /if -band images of thirty star-forming galaxies at z ~ 2 from the spectroscopic sample of 
iReddv et al.l (|2006T ) showing the mix of morphological types is comparable to the SMGs. As in Fig. [I] each image is 40 kpc at the redshift 
of the galaxy so that a direct comparison can me made. 
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Figure 4. Size versus redshift relation for SMGs (filled red circles) compared to the UV-SF (filled blue square) and high-z field (small 
black squares) comparison samples. The half light radii in the left hand panel are derived using HST /-band data whilst the right hand 
panel are derived from the NICMOS if-band (F160W) imaging. The dashed lines show that the medians of the distribution and illustrate 
that the SMGs and comparison samples have comparable half light radii in both the /- and if-bands. On the right hand side of each 
sub-panel we show the cumulative histograms of each distribution which shows that in the rest-frame UV the SMGs have marginally 
larger half light radii on average than the comparison samples, but essentially indistinguishable distributions in the rest-frame optical. 



exposure was corrected for a pedestal offset, and then mo- 
saiced together using the CALNICB task in IRAF. Unfortu- 
nately several exposures were affected by the South Atlantic 
Anomaly (SAA), and extra processing steps were requirecQ. 
The final images appear very fiat and have very low cosmic 
ray contamination. The observed / and i/-band images of 
each SMG in this sample are shown in Fig. [T] 

2.2 Comparison Sample 

To construct a comparison sample, we exploit the extensive 
spectroscopy in GOODS-N. To act as a high redshift "f ield" 
sample, we use ^2100 galaxies from lBarger et all (|2QQ8h be- 
tween z — 1-3.5 of which 330 also lie within the GOODS-N 
NICMOS imaging (Conselice et al. 2010 in prep). To en- 
sure a fair comparison, we restrict this high-redshift field 
sample to have the same redshift distribution as the SMGs 
(we note that the final high-redshift field sample has an I- 
band magnitude distribution of 1—22. 78=b 1.05; for reference, 
the SMGs have i=23.37±1.6). We hereafter refer to this 
sample as the "high-z field", but caution that the selection 
function for this sample is extremely complex since it com- 
prises a highly incomplete mix of UV/optical, mid-infrared 
and X-ray selected galaxies. We therefore also restrict the 
analysis specifically to a more homogeneously selected sam- 
ple of UV/optically selected galaxies (star- forming BX/BM 
and LBGs) with secure spectroscopic redshifts (in the range 
z -1.6-3.5) from lReddv et al.l (120061 ) (hereafter called UV- 
SF galaxies). This sample contains 228 galaxies which are 
in the GOODS-North ACS /-band imaging of which 53 are 
also have also been observed with NICMOS. We also re- 
move from the UV-SF sample one galaxy which has a bolo- 
metric luminosity greater than 1O 12 L0 and note that the 
median bolometri c luminosity of th e remaining sample is 
Uoi -1O 114 L (jReddv et all l2006h which is an order of 
magnitude lower than the typical bolometric luminosity of 

1 http: //www. stsci . edu/hst/nicmos/tools/post_SAA_tools .html 



the SMGs. In Fig.[2]and [3]we show true-colour HST images 
of the SMGs and UV-SF galaxies to demonstrate the visual 
mix of morphological types is comparable. 



2.3 Galaxy Sizes and Morphologies 

To quantify the galaxy morphologies (and measure their 
physical scale), we first calculate the Petrosian and half-light 
radii in the optical and near- infrared. The Petrosian radius 
is defined by ^^=1.5x^=0.2 where 77 = 0.2 is the radius 
(r) at which the surface brightness within a n annulus at r is 
one f i fth of the surface bright ness within r (|Conselice et al.1 
2000; Chapman et al . 2003b). This provides a measure of 
the size which does not depend on isophotes. The half-light 
radius, r^, is then defined as the radius at which the flux is 
one-half of that within r pe t- 

Given the small angular sizes of these objects, and their 
apparently complex morphologies, it is difficult to apply a 
standard morphological analysis along the lines of a "Hub- 
ble tuning fork" . Fortunately, statistics have been developed 
to help characterize high redshift galaxy m orphologies. We 
concentrate on the G ini and Asymmetries ([Abraham et al.l 
2003: IConselicl 120031 ) although we note that a number of 
other parameterisations of galaxy morphologies have been 
developed (such as concentration and dumpiness as part of 
the CAS system; IConselicd I2003T ). However, since the con- 
centration parameter represents the scale of the galaxy, and 
the dumpiness defines the spatial light distribution, to keep 
the analysis concise here we concentrate on the half-light 
radius and Gini co-efficient which are similar parameterisa- 
tions. 

The Gini coefficient is a statistical tool, originally de- 
veloped for economics, which determines the distribution of 
wealth within a population. Higher values indicate a very 
unequal distribution (G=l indicates all of the flux is in one 
pixel), whilst a lower value indicates it is more widely dis- 
tributed (G=0 suggests the flux is evenly distributed). The 
value of G is defined by the Lorentz curve of the galax- 
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Figure 5. Comparison between the structural parameters for 
SMGs compared to star-forming field galaxies. Top: Gini versus 
Asymmetry measured in the observed /-band for the SMGs show- 
ing that the median asymmetry for the SMGs and comparison 
samples is comparable, but that there is an offset between the 
Gini co-efficients of AC ~0.1 between SMGs and the comparison 
samples. The dashed line illustrates the correlation between Gini 
and Asymmetry in local galaxies (E-Sd morphological types) from 
Lotz et al. (2004). The cumulative histograms on each axis show 
that there is a significant difference in the rest-frame UV Gini 
coefficient for the SMGs compared to the comparison samples, 
but the asymmetries are comparable to the comparison samples. 
Centre: if-band Gini versus Asymmetry for the SMGs and field 
star-forming galaxies. As in the observed /-band the asymmetry 
is indistinguishable between the populations. However, the me- 
dian Gini coefficient for the SMGs in only ~0.03 larger than the 
field population. The dashed line denotes the same correlation as 
shown in the top panel. The cumulative histograms show that 
there is a subtle difference between the Gini co-efficient in the 
rest-frame optical for the SMGs and comparison samples, but the 
asymmetries are indistinguishable. Bottom: AG versus AA for 
the SMGs and field galaxies showing that the both the SMGs 
and field populations tend to prefer larger AG values. 



ies light distribution, which does not take into considera- 
tion any of the spatial information. The value G is derived 
by first sorting the pixel values within the Petrosian aper- 
ture, and then summed over a cumulative distribution (see 
lAbraham 6^111120031 ). 

Since the Gini coefficient removes all spatial informa- 
tion, to study the structures and morphologies of our sample, 
we also compute the galaxy asymmetry (A). The asymmetry 
of a galaxy is measured by minimising the subtraction of a 
galaxy image from itself after rotating the original image by 
180°. The minimisation takes in to account the uncert ainty 
in deriving the galaxy center fsee lConselice et all (|2QQ8h and 
references therein for a detailed discussion). Briefly, a lower 
value of A suggests that the galaxy is more symmetrical (e.g. 
elliptical galaxies), whilst higher values of A suggest highly 
asymmetric systems, usually found in spiral galaxies or ma- 
jor mergers (for reference, in the rest-frame UV/optical, an 
asymmetry value in excess of A=0.30 has been found to be a 
robust indication of a major merger at z=0 IConselice et al.1 
l2003aT ) 

To derive the asymmetry for the galaxies in our sam- 
ple, we first extract a 10" x 10" thumbnail of a galaxy image 
from the background subtracted frame. We mask all objects 
except for the galaxy image using the sextractor segmenta- 
tion map, replacing these regions with values taken from the 
correct noise properties as measured from the sky. We then 
measure the petrosian radius, half light radius and Gini co- 
efficient (r V et ,fh & G). To derive the asymmetry, we first 
derive a crude centre for the galaxy by finding the luminos- 
ity weighted centroid of the galaxy light distribution. We 
extract, rotate and subtract the galaxy image about this 
center, but then iterate and allow the center to vary over 
the entire half-light radius. This function is minimised to 
provide the integer pixel centre of the galaxy. Operationally, 
after the initial centre is found, the asymmetry is computed 
again for centres at the surrounding eight points in a 3x3 
grid. We use a distance of 0.1 pixels, corresponding to ap- 
proximately 0.1% of the galaxy half light radius and use bi- 
linear interpolation on the shifted image. If the asymmetry 
parameter at the centre is lower than any of the surround- 
ing pixels, then the asymmetry parameter is taken as this 
value. If the central pixel does not give the asymmetry min- 
imum, then the procedure is repeated with the new centre 
where the minimum was found. This process repeats until 
the minimum is found. 



3 ANALYSIS AND RESULTS 
3.1 Galaxy Sizes 

First, we compare the half light sizes of the SMGs and com- 
parison samples. As Fig. [4] shows, the SMGs have an opti- 
cal median half light radii of r opt (SMG)=2.3±0.3 kpc, which 
significantly overlaps with the half-light radii of UV-SF 
galaxies (r° h pt (UV - SF) = 1.9±0.2 kpc), and the field popu- 
lation (r^ pt (high — z field)=2.0=b0.1 kpc) (here we quote the 
error on the mean and we note that in all cases the la 
scatter in the distribution is is 1.3-1.7 kpc). These rest- 
fram e UV SMG ha l f light radii are similar to those found 
bv IChapman et al.1 (|2003bT ) who derive r^=2.8±0.3 kpc for 
a small sample of SMGs. 
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Turning to the observed if -band, the SMGs have a 
median half light radii of r% ir (SMG)=2.8±0.4 kpc which 
is in good agreement with both the high-z field popula- 
tion, and the UV-SF sample which have r^ r =2.6±0.2 and 
2.5±0.2kpc respectively. Thus it appears that the sizes of 
the SMGs are not systematically larger than the lower- 
luminosity star-forming galaxies at the same epoch in ei- 
ther f- or if -bands (Table 2), but the near-infrared half 
light radii tend to be systematically larger (~0.5 kpc) in all 
three galaxy populations. Finally, in Fig. [4] we show the op- 
tical and near-infrared sizes for all three galaxy samples as 
a function of redshift from z =0.5-3.5. Binning the samples 
into Az=0.5 bins, none of the three samples show systematic 
trends of half light radius with redshift. 

3.2 Gini and Asymmetries 

In Fig. [5] we show the I- and if-band Gini and 
Asymmetry coefficients for the SMGs compared to 
the UV-SF sample and the high-z field population. 
In the observed i-band the SMGs have a median 
Gini value of G opt (SMG)=0.69±0.03. In contrast, the 
comparison samples have G opt (UV - SF)=0.59±0.02 and 
G opt (high - zfield)=0.60±0.05. Thus the SMGs have a sys- 
tematically larger optical Gini coefficient (AC ~0.1) than 
the comparison samples. In the observed if-band, the 
SMGs also have a slightly larger Gini co-efficient with 
G n * r (SMG)=0.56±0.02 which is AG = 0.03 larger than the 
UV-SF population. The differences in the Gini co-efficient 
between SMGs and UV-SF galaxies (particularly in the 
rest-frame UV) suggests that dustier galaxies have star- 
formation which is less uniform, and/or that they suffer from 
more structured dust obscuration. 

Turning to the asymmetries, in both the rest-frame UV 
and optical, the SMGs have comparable asymmetries as the 
high-^ field and UV-SF samples (A opt (SMG) =0.27±0.03 
and A nir (SMG) =0.25±0.02). In the asymmetries alone, it 
therefore appears that the SMGs are no more likely to ap- 
pear as major mergers in the rest-frame UV/optical than 
more typical (lower bolometric luminosity), high-redshift 
galaxies. 

3.3 Light Profile 

Given the high resolution and reasonable signal-to-noise of 
our data, we also model the surface brightness distribu- 
tion of the SMGs and comparison samples. We use GALFIT 
(|Peng et al.ll200 2h to model both the I- and if-band surface 
brightness distributions for the sample. To ensure robust 
results, we restrict the analysis t o galaxies brighter than 
Iab <23.8 and for H A b <24.5 fsee lCimatti et alJl2008h . For 
each galaxy image, the PSF which GALFIT uses during the 
fitting process was generated with tinytim. In all cases, we 
used the Sersic profile to model the galaxy surface bright- 
ness profile, allowing the Sersic index, n, to vary between 
n=0.5 and n=6. For n=l and n=4 the Sersic profile reduces 
to an exponential or de Vaucoulers profile respectively. Disk 
dominated galaxies typically have low Sersic indices (n <2), 
whilst bulge dominated galaxies tend to have higher Sersic 
values (e.g. n >2). 

Individually, the SMGs have a wide range of Ser- 
sic indices in both the f- and if-bands, ranging from 



n ^0.5 to n=4.5, with a median m(SMG) = 1.8±1.0 and 
7iij(SMG) = 1.4±0.8 (the errors denote the error on the 
mean; the errors on individual measurements are typically 
±1.0). Similarly, the UV/optical star- forming comparison 
sample has a median Sersic index of m (UV-SF) = 1.2±0.8 
and riH (UV-SF) = 1.5±0. 7, although also with a comparably 
wide spread as the SMGs. 

An alternative test of the average light profile can be 
made by stacking the individual images. To achieve this, 
we first normalise the light distribution using the total 
flux within the Petrosian aperture (we only include galax- 
ies brighter than I<23.8, if <24.5 as above), and use the 
centre as defined by the Asymmetry minimisation proce- 
dure. In Fig. [6] we show the stacked two-dimensional im- 
ages, as well as the average light profiles for the SMGs and 
comparison samples. We also show the best-fit model and 
residuals after subtraction of the best fit profile. Using gal- 
fit to model the light profile (accounting for the PSF us- 
ing tinytim) the resulting best-fit model in the f-band has 
nj (SMG)=2.6±0. 5, whilst the if-band light profile is best 
described with n H (SMG)=2.0±0.5. Although the SMGs are 
extended on scales greater than the PSF in both f- and if- 
band (typically r^ >0.25 // ;2 kpc), it is likely that the spatial 
light profiles are sensitive to the PSF. We therefore con- 
volve the f-band image with the NICMOS PSF and the 
if-band image with the ACS PSF and re-fit the f- and if- 
band images with the convolved PSF respectively, obtaining 
nj(SMG)=2.5±1.0 and n H (SMG)=3.0±0.9. Since this extra 
smoothing produces Sersic indices within la of the previous 
measurements, we conclude that the Sersic indices give a rea- 
sonable estimate of the true Sersic indices of the galaxies. 
We also note that removing those SMGs with possible AGN 
contributions (Table 1), the Sersic index increases in both 
cases by An=0.2±0.3 (where the error accounts for the lower 
number of galaxies in the stack). Thus, the best fit SMGs in 
both I- and if-bands is n ~ 2.0-2.5 which is consistent with 
bulge-dominated galaxies. The same stacking procedure for 
the UV-SF comparison sample yields a slightly smaller Sersic 
index with m(UV-SF) = 1.8±0.3 and n H (UV-SF) = 1.5±0.3, 
significantly overlapping the Sersic values in both bands for 
the SMGs. 



4 DISCUSSION 

Using deep HST I and if-band imaging, we have performed 
a quantitative morphological analysis of a sample of 25 spec- 
troscopically confirmed sub-mm selected galaxies. We mea- 
sure the sizes of the SMGs, and derive typical half light radii 
ofr° p *(SMG)=2.3±0.3kpc and r£ ir (SMG)=2.8±0.4 kpc. We 
find that the SMGs have comparable sizes to UV-SF galax- 
ies, and the general z = 1-3 field population in both bands. 
The sizes we derive are slightly larger, but comparable to the 
gas sizes measured using IRAM/PdBI millim eter interferom- 
etry fn /o=1.8±0.8 kpc: [Tacconi et al.ll2008h . The similarity 
between the CO and observed if-band (stellar) sizes suggest 
that the stars and dense gas are most likely co- located. 

Although the measured sizes of SMGs and UV-SF 
samples are comparable, previous results have shown that 
there are significant differences in stellar and dynam- 
ical masses, with SMGs having potential well depths 
much greater than optically selected star-forming galax- 
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Table 1. 

Log of SMG sample with NICMOS imaging 



ID 


Short Name 


z 


TvDe 


F160W 


r nir 
' pet 


nir 
1 h 


r opt 
' pet 


°vt 

' h 












V K P C J 






V K P C J 


SMM J030227.73+000653.5 1 


CFRS03-15 


1.408 


SB 


20.73 + 0.02 


14.4+2.0 


2.8+0.4 


13.1 + 1.5 


2.3+0.5 


SMM J105158. 02+571800.2 


LOCKMAN-03 


2.239 


SB 


21.58 + 0.04 


9.9+1.3 


3.7+0.3 


7.8+1.3 


3.3+0.4 


SMM J105200.22+572420.2 


LOCKMAN-08 


0.689 


SB 


22.98 + 0.08 


5.3+0.5 


2.1+0.4 


5.0+1.3 


2.4+0.4 


SMM J105230.73+572209.5 


LOCKMAN-06 


2.611 


SB 


22.11 ±0.05 


10.2+0.5 


4.1+0.4 


15.9+2.5 


8.2+0.9 


SMM J105238. 30+572435. 8 


LOCKMAN-02 


3.036 


SB 


22.43 + 0.08 


8.6+0.6 


3.2+0.4 


5.7+2.0 


3.7+1.0 


SMM J123553. 26+621337.7 


HDFN-082 


2.098 


SB 


24.30 + 0.21 


3.1+0.2 


1.6+0.4 


3.7+0.4 


1.9+0.5 


SMM J123600. 10+620253. 5 2 


HDFN-092 


2.701 


SB 


23.75 + 0.16 


5.0+1.5 


2.9+0.4 


3.6+1.2 


2.4+0.4 


SMM J123600. 15+621047.2 


HDFN-093 


1.994 


SB 


22.64 + 0.08 


4.2+2.0 


2.3+0.4 


5.0+0.8 


1.5+0.4 


SMM J123606. 85+621021. 4 


HDFN-105 


2.509 


SB 


21.61 + 0.03 


4.2+0.5 


1.7+0.5 


7.6+0.5 


2.7+0.4 


SMM J123616.15+621513.7 


HDFN-127 


2.578 


SB 


23.59 + 0.10 


9.7+2.0 


2.1+0.6 


7.6+1.2 


2.0+0.4 


SMM J123622.65+621629.7 


HDFN-143 


2.466 


SB 


23.32 + 0.09 


9.7+0.4 


3.9+0.4 


6.1+2.0 


2.3+1.0 


SMM J123629.13+621045.8 


HDFN-153 


1.013 


SB 


21.23 + 0.03 


12.0+0.9 


4.5+0.4 


13.8+2.5 


5.6+0.9 


SMM J123632.61+620800.1 


HDFN-161 


1.993 


AGN 


22.74 + 0.07 


5.0+0.4 


2.0+0.5 


2.5+0.5 


0.9+0.4 


SMM J123635. 59+621424.1 


HDFN-172 


2.005 


AGN 


21.59 ± 0.03 


9.2+0.4 


3.7+0.4 


10.6+1.2 


1.8+0.5 


SMM J123701. 59+621513. 9 


GN17 


1.260 


SB 


21.50 + 0.05 


12.5+0.4 


5.3+0.4 






SMM J131201. 17+424208.1 


SA13-332 


3.405 


AGN 


22.87 + 0.08 


3.3+0.4 


1.3+0.4 


2.2+0.4 


0.4+0.4 


SMM J131225.20+424344.5 1 


SA13-516 


1.038 


SB 


21.46 + 0.03 


10.3+1.2 


3.9+0.4 


8.9+0.5 


2.7+0.4 


SMM J131232. 31+423949. 5 


SA13-570 


2.320 


SB 


22.85 + 0.09 


7.5+0.9 


2.5+0.5 


3.6+1.5 


1.1+1.0 


SMM J141741.81+522823.0 1 


CFRS14-13 


1.150 


AGN 


18.67 + 0.02 


5.6+1.5 


< 1 


6.6+2.0 


1.0+0.4 


SMM J141800.40+522820.3 1 


CFRS14-3 


1.913 


SB 


22.74 + 0.07 


3.7+0.6 


< 1 


8.3+0.3 


2.3+0.4 


SMM J163631. 47+405546. 9 


ELAIS-13 


2.283 


AGN 


24.47 + 0.18 


8.6+0.6 


4.2+0.4 


7.3+2.0 


4.3+1.0 


SMM J163639.01+405635.9 


ELAIS-07 


1.495 


SB 


23.88 + 0.09 


8.3+0.7 


3.0+0.4 


12.7+3.0 


3.6+0.8 


SMM J163650.43+405734.5 


ELAIS-04 


2.378 


SB/ AGN 


21.85 + 0.04 


9.7+0.6 


3.4+0.4 


5.0+1.5 


1.5+0.4 


SMM J163658.78+405728.1 


ELAIS-08 


1.190 


SB 


21.20 + 0.03 


9.9+0.6 


3.5+0.4 


6.8+2.0 


1.9+0.4 


SMM J163704.34+410530.3 


ELAIS-01 


0.840 


SB 


22.36 + 0.03 


6.9+1.2 


2.3+0.5 


6.1+0.5 


2.4+0.5 



Table 1. Notes: The ID's are taken from Chapman et al. 2005. ID's marked with a 1 denote galaxies observed with WFPC2 in the /-band 
whilst the galaxy marked by 2 d enotes observations wi th ST IS. The starburs t (SB) versus AGN classification is taken from optical and 
near-infrared spectroscopy from Chapm an et al.l {2005) and Swinban k et al.l (|2004T ). r pe t and denote pertrosian and half light radii 
respectively. SMM J123635. 59+621424.1 (GN17) is optically faint and so we have not attempted to derive petrosian and half light radii. 



Table 2. 

Median Morphological Parameters 



Sample 


r pS( k P c ) 


r£&(kpc) 


r° p *(kpc) 


r™ r (kpc) 


Qopt 




p^opt 


pnir 


AG 


AA 


SMGs 
UV-SF 
High- z field 


6.9+0.7 
6.2+0.2 
7.7+0.2 


7.7+0.6 
6.6+0.3 
7.8+0.2 


2.3+0.3 
1.9+0.2 
2.0+0.1 


2.8+0.4 
2.6+0.2 
2.5+0.2 


0.69+0.03 
0.59+0.02 
0.60+0.05 


0.56+0.02 
0.53+0.02 
0.55+0.01 


0.27+0.03 
0.29+0.02 
0.24+0.10 


0.25+0.02 
0.24+0.02 
0.24+0.01 


0.13+0.04 
0.06+0.03 
0.05+0.05 


0.02+0.03 
0.05+0.03 
0.00+0.10 



ies (e . g. ISwinbank et alJl2004l: lErb et al.ll2006l; iGenzel et al] 
2006; Forster Schreib e7et~al] 120061 : ISwinbank et al.1 [2006; 



Law et al.ll2007ah . If the regions sampled by the Hoj and CO 



emission lines are the same as those seen in the rest-frame 
UV/optical imaging, then this suggests that the SMGs have 
surface matter densities up to an order of ma gnitude larger 
than typical BX/BM and LBGs (see also iTacconi et al.l 
I2008T ). 

Using the asymmetry parameter to gauge the impor- 
tance of major mergers, we find that the SMGs and UV- 
SF comparison samples all have comparable asymmetries 
in the optical and near- infrared, with A = 0.27 =b 0.03 and 
A = 0.25 =b 0.02 in the observed /- and i7-band respec- 
tively. This is somewhat larger than typically measured for 
spiral or elliptical galaxies in the local Universe (typically 
A < 0.05), but comparable to the asymmetry measured in 



local ULIRGs (A = 0.35±0.1: IConselice et aT1 l2003a). Over- 
all, (and surprisingly) this suggests that in the the rest-frame 
UV/optical morphologies, SMGs are as equally likely to ap- 
pear as major mergers than lower- luminosity (and hence 
more quiescent) high-redshift star- forming galaxies. How- 
ever, the Gini co-efficients of the SMGs are systematically 
larger than the UV-SF galaxies, suggesting less uniform star- 
formation in the rest-fram e UV, and poss ibly structured 
dust obscuration (see also lLaw et al.ll2007bh . 

The striking result that the SMGs have comparable 
sizes and asymmetries as other high-redshift populations is 
at odds with the traditional picture of SMGs which have 
shown that SMGs are ex tended starbursts which are a result 
of major mergers (e.g. iGreve et al. 2005; ISwinbank et al. 
120061 : ITacconi et al]|2008l ). However, many of these stud- 



ies have concentrated on other multi-wavelength data (such 
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Figure 6. Average spatial light distribution of the SMGs and UV-SF galaxies in our sample in the observed I- and if-bands in the left 
and right hand panel respectively. In both panels, the flux scale is arbitrarily normalised and the I- and H- band profiles are offset in 
flux scale for clarity. The solid curves in both represent the best fit profiles. Inset: the optical and near-infrared stacks of the galaxies 
with the best-fit model and residuals. 



as kinematic studies through CO or Ha). Indeed, recently 
Ivison et al. (2010) use deep EVLA radio imaging of the 
CO (1-0) to show that the archetypal lensed SMG behind 
Abell 1835 (L1L2) has a large, extended cold molecular 
gas reservoir which extends over 25kpc in projection. To 
reconcile this apparent contradiction, we measure the half 
light radius, Gini co-efficient and asymmetry of L1L2 us- 
ing the source-plane HST imaging, deriving half light radii 
of r^=3.0±0.3 and 2.8±0.4kpc, asymmetries of A=0.32 and 
0.30 and Gini co-efficients of G=0.82 and 0.88 in observed I- 
and i/-bands respectively (centered on LI in this system). 
Thus, the rest-frame UV/optical morphology of L1L2 ap- 
pears similar to the typical SMGs studies here, even though 
the kinematics and spatial extent of the gas in this system is 
complex. Our findings suggest that, whilst multi- wavelength 
studies (such as those carried out via radio, CO or Ha) 
suggest archetypal SMGs may mark the sites of complex, 
merging systems, this is poorly reflected in their rest-frame 
UV/optical morphologies (or, at least it is difficult to differ- 
entiate SMGs from other more quiescent high-redshift pop- 
ulations via rest-frame UV/optical morphologies alone). 

Finally, to test how the stellar densities of the SMGs 
compare to other high redshift populations we combine the 
size and stellar mass estimates and show the size versus 
stellar mass relation in Fig. [71 We use the latest estimates 
for the SMG stellar masses from Hainline et al. (2009a, b) 
(whi ch take account of the TP-AGB phase of stellar evolu- 
tion; iMarastonl Il998l ) and AGN contribution to the rest- 
frame near-infrared photometry. For the SMGs, the me- 
dian stellar masses, (M* = 1.4 ± 0.5 x 10 11 M ) suggest 
stellar densities which are a factor ^5x large r (on av- 
erage) than the UV-SF po pulation (see also iBorvs et al.l 
120051 : lAlexander et alJl2008h . but slightl y larger (by a fac - 
tor 2x) than local early type galaxies (She n et al 1 120031 ). 
We also compare the sizes and stellar den sities to other 
z ^1. 5 luminous r ed galaxies (LRGs) from ICimatti et al.l 
<|2008h : IZirm et al.1 (|2007h which (due to their moderately 
high stellar masses, ~1O 1O ' 5_11 ' O M , compact sizes rn — 
1.2 ± 0.2 kpc, and low space densities, ^10 -4 Mpc -3 ) have 
been proposed as direct SMG descendants. However, the 
stellar masses and sizes of the SMGs are, on average, both a 



factor ^2x larger than LRGs. Thus, unless the stars forming 
in the gas reservoir ultimately have a very different spatial 
distribution than the gas itself (which is unlikely given the 
similarity between the observed stellar and gas sizes), it is 
difficult to reconcile the small physical sizes of luminous red 
galaxies at z ~1.5 with an evolutionary sequence with SMGs 
which already appear to be a factor ~ 2x larger measured 
at the same rest-frame wavelength at z ~ 2.5. 



5 CONCLUSIONS 

We have undertaken the first large near-infrared morpholog- 
ical analysis of SMGs. We find that the SMGs have compa- 
rable sizes to UV-SF galaxies at the same epoch (BX/BM 
and LBGs) in both /- and i/-bands, and (surprisingly) 
with comparable asymmetries. However, we find that the 
SMGs have systematically larger Gini co-efficients (particu- 
larly in the observed /-band) than UV-SF galaxies at the 
same epoch suggesting less uniform, high intensity star- 
formation in the rest-frame UV, possibly reflecting struc- 
tured dust obscuration (see also lLaw et a l. 2007b). Overall 
our results suggest that SMGs are no more likely to appear 
as major mergers in the rest-frame UV/optical than more 
typical, lower- luminosity high-redshift galaxies (such as Ly- 
man break galaxies or BX/BMs). However, the differences 
in the Gini coefficients between populations suggests that 
the dustier SMGs have star- format ion which is less uniform 
(and/or that they suffer from more structured dust obscu- 
ration) . 

We also show that most of the SMGs have observed H- 
band light profiles which are better fit with that of a spheroid 
galaxy light distribution. Indeed, stacking the galaxies we 
find that in the observed if-band, an n~2 Sersic index pro- 
vides a better fit to the spatial light profile than an expo- 
nential disk model, suggesting that, whilst these galaxies are 
individually morphologically complex, the composite stellar 
structure of the SMGs reflects that of a spheroid/elliptical 
galaxy. However, we note that the same analysis of the UV- 
SF galaxies is statistically indistinguishable, with only a 
marginally lower Sersic index with n — 1.6 =b 0.3. 

The close similarity between the rest-frame UV and 
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Figure 7. The distribution of physical sizes (77J versus stellar mass for SMGs compared to K-band selected samples of luminous 
red galaxies at z ~l-2. Stellar masses are estimates using the iMarastonl (1998) stellar population libraries. The solid line shows the 
size-mass relation of early- type galaxies from SDSS by Shen et al.l (|2003[) (with d otted lines indicating the la scatter). The comparison 
sampl es comprise the GMASS passive, de nse galaxies from | Cimatti et all (|2008l ). passive galaxies selected from FIRES bv IZirm et al.l 
(2007), massive galaxies from MUNICS bv lTruiillo et al.l (2006). We also include the optical/UV selected star-forming galaxy samples by 
combining their sizes (this work) and median stellar mass ([Shaplev et al.ll2p05l). For clarity we spl it the SMGs into two samples: starburst 
(SB) and AGN where the classification is based on the spectroscopy of IChapman et al.l (2005); Swinbank et al ] (|2004h ; Takata et al. 
([20061 ). This figure shows that approximately half the SMGs have comparably high stellar mass densities as the passive, luminous red 
galaxies from GMASS, but with substantially larger sizes and stellar masses (by a factor ~2x in both cases). 



optical morphologies in the SMGs, and UV-SF galaxies 
suggests that both wavelengths are dominated by young, 
star-bursting compon e nts as well as dusty r egions (see also 
iDickinson et al.l l2000l : IPapovich et al.ll2005h . These results 
are in contrast to local studies of similarly luminous LIRGS 
and ULIRGs in the local Universe which have been shown 
to have very different Hubble types from the rest-fr ame UV 
and optical wavelengths (e.g. iGoldader et al.l 120020 . possi- 
bly suggesting fundamental differences between starbursts 
at z = and z ~ 2. 

Although the sizes and structural properties of the 
SMGs and UV-SF star-forming galaxies are similar, pre- 
vious work has shown that the dynamical masses of 
SMGs are up to an order o f mag n itude larger than UV- 
SF galaxies (e . g. lErb et all I2QQ3I : ISwinbank et al.l l2004|: 
Erb et alJlioOft iGenzel et al.|l2006l: iForster Schreiber et all 
20061 : ISwinbank et all l200d lLaw et all l2007ah .. This sug- 
gests that the intense star- format ion within SMGs does not 
represent an evolutionary sequence in which typical UV- 
SF galaxies undergo intense star- format ion (either through 
merging or through secular processes), but rather it is the 
availability of large gas reservoirs within already massive 
galaxies that allow the SMG phase. 

Finally, we also investigate the size-stellar mass re- 
lation of SMGs in order to test whether the SMGs may 
represent progenitors of the luminous red galaxies seen at 
z ~1.5. We combine estimates of the size and stellar masses 
to show that approximately half of the sample have stel- 
lar mass densities comparable to those derived for luminous 
red galaxies (eg. ICimatti et al.l l2008h . However, we show 
that the median size of the SMGs in the observed near- 
infrared (ru = 2.3 =b 0.3 kpc) is larger than that of the lu- 
minous red galaxies, which have a median half light radius 



of tk — 1-2 =b 0.2 kpc. We also show that the median stellar 
mass of the SMGs is also a factor ~ 2x larger, thus suggest- 
ing that the luminous red population at z ~ 1.5 are unlikely 
to be direct descendants of the SMG population unless the 
new stars formed in the SMG starburst ultimately have a 
very different spatial distribution from their gas reservoirs, 
which seems unlikely given the similarity between the CO 
and UV/optical sizes. 

Overall, our results suggest that rest-frame UV and op- 
tical morphologies of high-redshift galaxies are essentially 
decoupled from other observables (such as bolometric lumi- 
nosity, stellar or dynamical mass). Alternatively, the physi- 
cal processes occuring within the galaxies are too complex to 
be simply characterised by the rest-frame UV/optical mor- 
phologies. It may be that at significantly longer wavelength 
structural differences will appear, but this will have to wait 
for rest-frame near-infrared imaging with James Webb Space 
Telescope, (JWST). Alternatively, high resolution kinemat- 
ical studies on sub-kpc scales (e.g. with near-infrared IFUs 
which can now be carried out from the ground using adap- 
tive optics) may offer the most direct route to probing the 
differences between high-redshift galaxy population as the 
dynamics, distribution of star- format ion and metallicity gra- 
dients will reflect differences in the triggering mechanism 
and mode of star- format ion at high-redshift. 
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